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INTRODUCTION 


In  gun 8  with  a  sliding  breech  mechanism,  breech  ring  failures  have  been 
observed  originating  from  the  lower  fillet  in  the  vicinity  of  the  contact 
region.  The  observations  indicate  that  high  tensile  stress  produced  by  stress 
concentration  at  the  fillet  was  responsible  for  the  failure.  In  order  to 
reduce  the  chance  of  failure  and  extend  the  fatigue  life,  an  exploratory  study 
was  initiated  on  the  autofrettage  of  a  breech  mechanism.  The  technique  is 
based  on  the  production  of  beneficial  residual  stresses  through  coldworking  to 
counteract  the  high  operating  stresses  induced  by  firing. 

A  photoplastic  model  made  of  aluminum  block  and  polycarbonate  ring  was 
designed. 1  The  maximum  fillet  stresses  for  an  elastic  load,  as  well  as  an 
elastoplast ic  load,  were  determined  experimentally.  Residual  stress  resulting 
from  removing  maximum  test  load  was  calculated.  A  numerical  investigation  of 
the  photoplastic  model  was  made  by  using  NASTRAN. ^  The  numerical  results  are 
in  good  agreement  with  the  experimental  data  in  the  elastic  range  of  loading, 
but  not  in  the  plastic  range.  The  NASTRAN  program  stops  before  the  maximum 
test  load  is  reached,  thus  residual  stress  after  complete  unloading  cannot  be 
calculated. 

This  report  describes  a  numerical  technique  for  predicting  the  residual 
stresses  in  an  overloaded  breech  ring.  The  numerical  approach  used  is  the 
finite  element  method  based  on  the  incremental  stress-strain  matrix,  for  which 

l-Cheng,  Y.  F.,  "Photoplastic  Study  of  Residual  Stress  in  an  Overloaded  Breech 
Ring,"  Technical  Report  ARLCB-TR-78018,  November  1978. 

^Chen,  P.  C.  T.  and  O’Hara,  G.  P.,  "NASTRAN  Analysis  of  a  Photoplastic  Model 
for  Sliding  Breech  Mechanism,"  Proceedings  1979  Army  Numerical  Analysis  and 
Computer  Conference,  pp.  237-254. 
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a  computer  program  is  developed.  '  The  material  behavior  is  characterized  by 
the  von  Mises1  yield  criterion,  Prandtl-Reuss  flow  equations  and  isotropic 
hardening  rule.  A  piecewise  linear  representation  for  the  stress- strain  curve 
is  used.  The  numerical  results  of  the  stresses  in  all  elements  are  obtained 
for  four  contact  conditions.  The  location  and  magnitude  of  the  maximum  fillet 
stress  are  determined  as  a  function  of  loading.  Residual  stresses  resulting 
from  elastic  unloading  are  calculated.  A  satisfactory  agreement  between 
numerical  and  experimental  results  has  been  reached. 

METHOD  AND  PROGRAM 

The  finite  element  method  used  here  is  of  the  incremental  type  since 
solutions  to  problems  involving  nonlinear  material  are  best  obtained  by  solv¬ 
ing  a  sequence  of  linear  problems  associated  with  an  incremental  application 
of  loading.  The  formulation  of  the  governing  matrix  equation  is  developed  by 
using  the  incremental  stress-strain  matrix  in  conjunction  with  the  stationary 
principle  for  elastic-plastic  solids.  Following  the  procedure  outlined  in 
reference  3,  we  have  developed  a  finite  element  computer  program  for  solving 
two-dimensional  elastic-plastic  problems  of  rotational  symmetry.  The  element 
type  is  either  triangular^  or  quadrilateral^,  and  both  were  implemented  in 

3Yamada,  Y., Yoshimura,  N. ,  and  Sakurai,  T.,  "Plastic  Stress-Strain  Matrix  and 
Its  Application  for  the  Solution  of  Elastic-Plastic  Problems  by  the  Finite 
Element  Method,"  International  Journal  of  Mechanical  Science,  Vol.  10,  pp. 
343-354,  1968. 

^Chen,  P.  C.  T.,  "Elastic-Plastic  Solution  of  a  Two-Dimensional  Tube  Problem 
by  the  Finite  Element  Method,"  Transactions  of  Nineteenth  Conference  of  Army 
Mathematicians,  ARO  Report  73-3,  pp.  763-784,  1973. 

3Chen,  P.  C.  T. ,  "Numerical  Solution  of  Gun  Tube  Problems  in  the  Elastic- 
Plastic  Range,"  Proceedings  1977  Army  Numerical  Analysis  and  Computer 
Conference,  pp.  423-439. 
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analyzing  gun  tube  problems.  For'  the  present  problem,  the  plane  stress  case 
was  implemented.  The  field  displacement  is  based  on  the  linear  displacement 
function  in  a  triangular  element.  Small  strain  condition  is  used  to  derive 
the  kinematical  matrix  [ B] .  The  material  behavior  is  characterized  by  the  von 
Mises'  yield  criterion,  Prandtl-Reuss  flow  equations  and  isotropic  hardening 
rule.  A  piecewise  linear  representation  for  the  effective  stress-strain  curve 
is  used.  As  derived  in  reference  3,  the  incremental  form  of  the  constitutive 
matrix  [ D]  for  an  elastic-plastic  element  (a  >  aQ)  can  be  written  as 

{Aa}T  =  [D]  {Ae} 

{Ao}t  -  Aax,  Aay,  ATxy 
{ Ae}t  -  Aex,  Aey,  Ayxy 

and 


[D] 


a3 


Oy’^A! 


-ax,°y,+2vAi  ,  ax'2+2A1 


SYM. 


Ox’+voy' 

1+v 


^xy  » 


Oy'+VOjj' 


A2  2H'(l-v)a 
1+v  xy  *  2(l+v)  9E 


(1) 


where  E,  v  are  Young's  modulus,  Poisson's  ratio,  respectively. 


^Yamada,  Y.,  Yoshimura,  N. ,  and  Sakurai,  T. ,  "Plastic  Stress-Strain  Matrix  and 
Its  Application  for  the  Solution  of  Elastic-Plastic  Problems  by  the  Finite 
Element  Method,"  International  Journal  of  Mechanical  Science,  Vol.  10, 
pp.  343-354,  1968. 
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2H'  -  Txy  „  „ 

Ai  - - a2  + - ,  Ao  =  ax' 2  +  2va'av'  +  av' 2 

1  9E  1+v  Z  y  y 


A3  =  2(1-v2)A1  +  A2  ,  ax  -  (2ax-ay)/3  , 

-  -p 

Oy'  =  (2ay-ax)/3  ,  H'  =*  Aa/Ae  , 


a  =  [ox2  +  Oy2  -  axay  +  3TXy2]l/2 
AeP  =  2/3  [(AePx)2  +  (Ae  P)2  +  (AePz)2  +  -  (Ay^)2]1/2 


(2) 


and  0 q  is  the  yield  stress  in  simple  tension  or  compression.  For  an  elastic 
element  in  which  a  <  aQ,  the  constitutive  matrix  is 


[D] 


E 


1-v^ 


1  v  0 
v  1  0 

0  0  (l-v)/2 


(3) 


The  stiffness  matrix  [k]  for  a  triangular  element  is  represented  by 

[k]  -  [B]T[D][B]tA  (4) 

where  t  denotes  the  thickness  and  A  is  the  area  of  the  triangular.  Four  tri¬ 
angular  elements  are  used  to  construct  a  quadrilateral  element  as  shown  in 
reference  5.  The  element  stiffness  matrix  is  determined  by  static  condensa¬ 
tion  and  the  kinematical  matrix  is  formed  by  averaging  at  the  centroid  of  each 
element . 


^Chen,  P.  C.  T. ,  "Numerical  Solution  of  Gun  Tube  Problems  in  the  Elastic- 
Plastic  Range,"  Proceedings  1977  Army  Numerical  Analysis  and  Computer 
Conference,  pp.  423-439. 
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The  computer  program  developed  for  the  plane  stress  elastic-plastic  prob¬ 
lems  was  a  modification  of  the  previous  program  for  the  axisymmetric  case.^>^ 
The  overlay  feature  was  utilized  for  reducing  the  core  storage  requirement*  A 
disk  was  used  as  an  auxiliary  storage  device  for  the  current  element  informa¬ 
tion  and  for  solving  the  large  system  of  equation  by  partition*  A  tape  may  be 
used  to  store  the  final  results  for  output  plotting  and  also  for  restarting 
this  program  from  a  point  of  completion  of  a  given  loading  sequence*  The  com¬ 
puter  used  was  an  IBM  360  Model  44. 

MODEL  AND  LOADING 

A  two-dimensional  photoplastic  model  of  aluminum  block  and  polycarbonate 
ring  was  designed. ^  The  breech  ring  was  made  of  0.12  inch  thick  LEXAN  plate 
and  the  top  of  the  ring  was  fixed.  This  material  has  a  Poisson’s  ratio  of 
0.38  in  the  elastic  state  and  a  limiting  value  of  0.5  in  the  plastic  state. 

The  stress-strain  (a-e)  curve  for  LEXAN  beyond  the  proportional  limit  can  be 
described  by  the  modified  Ramberg-Osgood  equation  in  the  following  form 

Ee/aB  -  a/aB  +  (3/7)(a/aB)n  for  aA  <  a  <  ac  (5) 

and  the  values  of  five  parameters  are  E  -  325  Ksi,  n  -  11.5,  aA  -  6.2  Ksi, 
aB  =»  8.7  Ksi,  Oq  =  9.576  Ksi,  where  E  is  Young’s  modulus,  n  is  a  parameter, 


l-Cheng,  Y.  F.,  "Photoplastic  Study  of  Residual  Stress  in  an  Overloaded  Breech 
Ring,”  Technical  Report  ARLCB-TR-78018,  November  1978. 

^Chen,  P.  C.  T. ,  "Elastic-Plastic  Solution  of  a  Two-Dimensional  Tube  Problem 
by  the  Finite  Element  Method,"  Transactions  of  Nineteenth  Conference  of  Army 
Mathematicians,  ARO  Report  73-3,  pp.  763-784,  1973. 

^Chen,  P.  C.  T. ,  "Numerical  Solution  of  Gun  Tube  Problems  in  the  Elastic- 
Plastic  Range,"  Proceedings  1977  Army  Numerical  Analysis  and  Computer 
Conference,  pp.  423-439. 
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is  the  proportional  limit,  Og  is  the  secant  yield  strength  having  a  slope 
equal  to  0.7  E  and  Oq  is  the  flow  stress  at  which  the  slope  of  the  stress- 
strain  curve  is  zero.  A  piecewise  linear  representation  for  the  effective 
stress-strain  curve  is  used  in  this  numerical  study  and  the  values  for  11  data 
points  other  than  the  origin  are 

(eif  a±  in  Ksi)  =■  (.018621,  6.),  (.020401,  6.5), 


■c 

(.022480, 

7.0), 

(.025158, 

7.5), 

(.028988, 

8.0), 

(.034934, 

8.5), 

(.044635, 

9.0), 

(.051679, 

9.25), 

(.060781, 

9.50) 

,  (.064043 

,  9.596),  (.1,  9 

.576)  . 

(6) 

The  dimensionless  slope  is  defined  as 

tOi  «  E”1  (oi-j-i-oi)/(ci+1-c^)  (7) 

and  we  can  calculate  (Hf/E)i  ■  ^i/(l-^i)  for  i  *  1,...,10. 

A  finite  element  representation  for  one  half  of  the  breech  ring  is  shown 
in  Figure  1.  The  other  half  is  not  needed  because  of  symmetry.  There  are  224 
grid  points  and  189  quadilateral  elements  in  this  model.  The  grids  1  thru  8 
are  constrained  in  x-direction  only  while  grids  217  thru  224  are  held  fixed. 
The  top  portion  of  the  breech  ring  is  omitted  because  it  is  believed  to  have 
little  effect  on  the  maximum  stress  information  near  the  lower  fillet.  This 
belief  has  been  confirmed  by  obtaining  the  elastic  solution  for  another  finite 
element  model  with  an  additional  70  quadilateral  elements  in  the  top  portion. 
The  difference  between  these  two  models  for  the  maximum  tensile  stress  is  only 
1.3  percent.  The  aluminum  block  is  regarded  as  rigid  and  the  load  is 
transmitted  to  the  ring  through  contact.  Initially,  the  block  is  in  full 
contact  with  the  ring.  As  the  load  increases,  a  gap  develops  in  the  central 
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portion.  The  width  of  the  central  gap  under  the  full  test  load  (2F  ■  1144 
pounds)  is  observed  experimentally  to  be  about  five  inches.  Our  elastic- 
plastic  finite  element  program  in  its  present  form  cannot  be  used  to  determine 
the  width  of  contact  and  the  force  distribution  as  functions  of  loading. 

Guided  by  the  experimental  information  on  the  width  of  central  gap,  we 
have  chosen  four  contact  conditions  in  this  numerical  investigation.  The 
points  of  contact  are  at  nodes  (33,  41,  49,  57)  for  case  1,  at  nodes  (41,  49, 
57)  for  case  2,  at  nodes  (49,  57,  65)  for  case  3  and  at  nodes  (57,  65)  for 
case  4.  The  width  of  contact  and  the  force  distribution  in  each  case  are 
assumed  to  remain  unchanged  during  loading  and  only  the  total  force  (2F)  is 
allowed  to  increase.  The  force  distribution  may  be  uniform  or  non-uniform  and 
both  types  have  been  used  in  this  study.  Initially  we  apply  a  small  force  to 
obtain  elastic  solution  and  the  total  force  (2F*)  required  to  cause  incipient 
plastic  deformation  is  calculated  by  using  the  von  Mises1  yield  criterion. 

Then  we  apply  the  additional  force  in  increments  until  the  maximum  test  load 
is  reached. 

The  load  increments  we  chose  are  non-uniform  because  our  experience  indi¬ 
cates  that  smaller  increments  should  be  used  as  plastic  deformation  becomes 
larger.  Our  choice  for  the  present  investigation  is  an  arithmetic  decreasing 
sequence  generated  automatically  in  the  program.  If  ten  incremental  steps  are 
used,  then 

1 

AFi  -  -  19, 18,... ,11, 10  (5.72-F*)  for  j  =  1,...,10  (8) 

J  145 

It  is  important  to  choose  a  proper  set  of  increments  in  order  to  obtain  good 
results  at  a  reasonable  cost.  To  reach  the  maximum  test  load  of  2F  -  1144 
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pounds,  we  used  11  increments  for  the  four  contact  conditions  and  22 
increments  for  the  second  and  third  cases.  The  difference  between  11  and  22 
incremental  loadings  for  the  maximum  tensile  stress  is  found  to  be  within  one 
percent.  The  restart  feature  of  the  program  has  been  used  for  the  second 
contact  condition  to  increase  the  total  force  from  1144  pounds  to  1300  pounds 
in  seven  additional  steps.  The  increments  are  again  non-uniform. 

RESULTS  AND  DISCUSSIONS 

The  numerical  results  of  the  stresses  in  all  elements  were  obtained  for 
the  overloaded  breech  ring  under  different  contact  conditions.  Some  of  them 
are  presented  below  in  Figures  2  thru  7.  The  major  principal  stresses  in 
elements  along  the  contact  region  and  fillet  are  shown  in  Figures  2  thru  5 
for  the  four  contact  conditions  under  uniform  load  distributions.  In  each  of 
these  figures,  we  presented  two  sets  of  data  corresponding  to  the  load  levels 
at  the  initial  yielding  (F*)  and  at  the  maximum  test  load  (F  -  572  pounds). 

The  incipient  plastic  deformation  first  occurs  in  element  99  for  the  first  two 
cases  and  in  element  50  for  the  last  two  cases.  The  values  of  load  level  for 
the  four  cases  are  F*  =  312.8,  321.9,  221.7,  169.2  pounds.  The  initial 
yielding  is  tensile  in  element  99  and  compressive  in  element  50.  It  is 
interesting  to  find  out  that  the  location  of  the  maximum  tensile  stress  is  in 
element  99  for  all  contact  conditions  under  either  uniform  or  non-uniform  load 
distributions  and  this  location  remains  unchanged  as  the  total  force 
increases.  Since  this  location  is  of  sufficient  distance  from  the  contact 
regions,  it  seems  to  suggest  that  Saint  Vennant's  principle  can  be  applied  to 
this  problem  in  the  elastic  as  well  as  plastic  range  of  loading. 
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The  values  of  the  maximum  tensile  stress  based  on  the  four  contact 


conditions  with  uniform  load  distribution  are  2029,  1972,  1939,  1936,  psi  at  F 
*  100  pounds  and  9532,  9361,  9174,  9119  psi  at  F  *  572  pounds.  It  should  be 
noted  that  the  stresses  in  the  finite  element  program  were  calculated  at  the 
centroid  of  each  element  but  only  one  principal  stress  at  the  boundary  was 
measured.  The  values  of  the  maximum  tensile  stress  at  the  fillet  based  on 
experimental  approach  are  2222  psi  at  F  *  100  pounds  and  9300  psi  at  F  ■  572 
pounds. 

For  the  purpose  of  comparison,  the  boundary  stress  is  determined  by 
extrapolation  using  the  calculated  results  for  those  elements  along  the  radial 
direction  through  element  99.  This  is  illustrated  in  Figure  6  for  the  second 
case  of  contact  condition.  Simmilar  figures  for  the  other  three  cases  are  not 
shown.  Three  curves  are  plotted  in  Figure  6  and  they  represent  the  major 
principal  stresses  for  three  load  levels  at  the  initial  yielding,  maximum  test 
load  and  complete  unloading  after  reaching  the  maximum  load.  The  residual 
stresses  after  complete  unloading  are  determined  by  assuming  that  the. 
unloading  process  is  purely  elastic.  Our  numerical  results  reveal  no  reverse 
yielding.  In  extrapolating  the  boundary  stress,  we  shall  remember  that  the 
maximum  tensile  stress  shall  not  exceed  the  flow  stress  of  9576  psi.  As  seen 
in  Figure  6,  a  comparison  between  the  numerical  and  the  experimental  results 
for  the  maximum  tensile  stress  indicates  that  a  satisfactory  agreement  has 
been  reached. 

In  Figure  7,  two  principal  stresses  as  well  as  the  residual  stresses  in 
element  99  are  shown  as  functions  of  loading.  Only  the  second  case  of  contact 
conditions  is  presented  here  for  illustration.  The  residual  stress  is 
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determined  by  assuming  a  purely  elastic  unloading  resulting  from  various 
stages  of  loading*  The  minor  principal  stress  (02)  is  found  to  be  a  nearly 
linear  function  of  loading  and  its  residual  value  is  very  small.  The  value  of 
the  principal  stress  angle  for  all  contact  conditions  and  for  all  load  levels 
is  found  to  lie  within  -26°  to  -27°  with  respect  to  the  x-axis.  As  can  be 
seen  in  Figure  7,  the  major  principal  stress  (o^)  and  its  residual  value 
increase  in  magnitude  as  the  total  contact  force  increases,  but  they  are  of 
opposite  signs*  Therefore,  the  production  of  compressive  residual  stress  in 
an  overloaded  breech  ring  is  beneficial  because  it  can  counteract  the  high 
operating  tensile  stress  induced  by  firing.  As  a  result  of  prestressing  dur¬ 
ing  the  process  of  manufacturing,  the  elastic  carrying  capacity  and  fatigue 
life  can  be  increased ,  and  thus  the  chances  of  breech  failures  can  be  reduced. 
Furthermore,  the  finite  element  program  developed  here  is  a  useful  tool  in 
redesigning  breech  rings  or  blocks. 
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Figure  1.  A  Finite  Element  Model  for  Breech  Ring 
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Figure  2.  Major  Principal  Stress  Along  the  Boundary- 
Elements  for  Case  1 
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Figure  3.  Major  Prlnoipal  Stress  Along  the  Boundary 
Elements  for  Case  2 
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Figure  4.  Major  Principal  Stress  Along  the  Boundary 
Elements  for  Case  3 
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Figure  5.  Major  Principal  Stress  Along  the  Boundary 
Elements  for  Case  4 
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Figure  6.  Determination  of  the  Maximum  Tensile  Stress  for  Case 
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Figure  7.  Stresses  in  Eleaent  99  as  functions  of  Contact  Force  for  Case 
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